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ABSTRACT 

We numerically model fragmentation of a gravitationally unstable gaseous disc un- 
der conditions that may be appropriate for the formation of the young massive stars 
observed in the central parsec of our Galaxy. In this study, we adopt a simple prescrip- 
tion with a locally constant cooling time. We find that, for cooling times just short 
enough to induce disc fragmentation, stars form with a top-heavy Initial Mass Func- 
tion (IMF), as observed in the Galactic Centre (GC). For shorter cooling times, the 
disc fragments much more vigorously, leading to lower average stellar masses. Thermal 
feedback associated with gas accretion onto protostars slows down disc fragmentation, 
as predicted by some analytical models. We also simulate the fragmentation of a gas 
stream on an eccentric orbit in a combined Sgr A* plus stellar cusp gravitational po- 
tential. The stream precesses, self-collides and forms stars with a top-heavy IMF. None 
of our models produces large enough co-moving groups of stars that could account for 
the observed "mini star cluster" IRS13E in the GC. In all of the gravitationally un- 
stable disc models that we explored, star formation takes place too fast to allow any 
gas accretion onto the central super-massive black hole. While this can help to explain 
the quiescence of 'failed AGN' such as Sgr A*, it poses a challenge for understanding 
the high gas accretion rates infered for many quasars. 
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1 INTRODUCTION 

There is no detailed understanding of how super-massive 
black holes (SMBHs) gain their mass, except that 
it must be mainly through gaseous disc accretion 
(|Yu k, Tremaine, 20021). The sta ndard accretion disc 
model (Shakur a Sunvaev, 1973l l shows that accre- 
tion discs must be quite massive and cold at "large" 
(sub-parsec in this context) distances from the SMBH. 
Many theorists pointed out that these discs should 
be unstable to self-gravity, and thus must form stars 
or giant planets by gravit ational f ragmentation (e.g., 



demonstrated that star formation cannot be quenched by 
stellar feedback, unless one is prepared to grossly violate 
constraints that we have from AGN spectra. 

i Goodman (2003 1) suggested that the feeding of SMBHs 
proceeds via direct accretion of low angular momentum gas 
that settles in a small scale disc. Such a disc would be 
too hot to allow star formation. The size of this "no-star- 
formation region" is R ^ 0.03 pc, and is weakly dependent 
on the SMBH mass. However, one realistically expects that 
an even larger amount of gas settles at larger radii where 
it could collapse and form stars, or it could accrete on the 



iPaczvnski. 197SKo l vkhalov Sunvaev. 198CShlosman Begelma^M^£g^Mm%^lite.f^^gjfejjlf^^ 200lh . 



This creates a dilemma for the field, as star forma- 
tion might be a dynamical, very fast process which 
may result in a complete transformation of the gas 
into stars. The SMBH would then b e starved of fuel. 
iGoodman f200dSirko Goodman (20031 ) have recently 
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still very much interesting. 

Current observations provide new im- 
petus for theoretical work on the topic. 
IWehner Harris (200dFerrarese et al. (20061 ) have shown 
that many galaxies host "extremely compact nuclei" 
(ECN) - star clusters located in the very central ^ 5 
parsec of galaxies. Closer to home, in the centre of our 
Galaxy, about a hundred massive young stars, mostly 

arranged into two stellar discs of ^ 0.05-0.5 parsec scale 

(|Levin k Beloborodov. 200dGenzel et al.. 2003lPaumard et al.. 2006l l. 
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circle the super-massive black hole named Sgr A*. These 
stars are suspected to have been formed in situ, as ev- 
idenced by the lack of stars outside the ^ 0.5 parse c 
region (^Paumard et al. . 200dNavakshin Sunvaev. 2005l l. 
Several observational facts are consistent with the hy- 
pothesis that these st ars formed in a massive gaseous disc 
(|Paumard et al.. 20061 ). Had this star formation event not 
happened, Sgr A* could have been accreting gas from the 
gaseous disc even now. Sgr A* thus failed to realise itself 
as an AGN because of the loss of gas to star formation 
([Navakshin Cuadra. 20051 ). 

A number of important gaps in our understanding 
of young stars near Sgr A* remain. The counter clock- 
wise disc appears to host stars on more elliptical orbits, 
and it is also geometrically thicker. The same disc also 
contains a puzzling "mini star cluster", IRS13E, that 
consists of more than a dozen stars and may be bound 
by an Intermed iate Mass Black Hole (IMBH) of mas s 
Mbh ^ 10^ M© (Pa umard et al.. 2 005Scho del et al.. 20051 ). 
The IMF of the observed stars must be top- 
heavy according to several lines of evidence 
([Navakshin Sunvaev. 2005lNa vakshi n et al.. 200dPaumard et 
which is not expected in the most basic model of a frag- 
menting disc, as the Jeans mass there is significantly 
sub-solar (e.g.. Levin 2006, Nayakshin 2006, but see also 
Larson 2006). 

In this paper, we discuss numerical simulations of star 
formation occurring in a gaseous disc around Sgr A*. Given 
the numerical challenges in the problem, we foresee that a 
reliable modelling of all the questions raised by the young 
stars near Sgr A* will require an extended effort of con- 
stantly increasing complexity. In this study we present nu- 
merical experiments with a locally constant cooling time. 
This allows a convenient comparison with previous analyti- 
cal and numerical works that predicted the conditions when 
fragmentation should take place. It might also form a basis 
for comparison with future work. 

Within our formalism with a locally constant cooling 
time, we find that (i) circular and eccentric discs alike can 
gravitationally fragment and form stars; (ii) the IMF of 
formed stars is a strong function of cooling time, becoming 
top-heavy for marginally star-forming discs; (iii) star forma- 
tion feedback is indeed able to slow down disc fragmenta- 
tion, as suggested by several earlier analytical papers, but it 
is not yet clear if it can alleviate the fueling problem of the 
SMBHs; (iv) our simulations do form some tightly bound 
binary stars but more populous systems (a la IRS13E) do 
not survive long. 

This paper is structured as follows. In Section [2l we 
discuss our simulation methodology and the basic conditions 
for disc fragmentation. Section [3] explains our sink particle 
approach to treat star formation in more detail. We then 
analyse the evolution after disc fragmentation and the IMF 
of the formed stars in Section U The sensitivity of our results 
to numerical resolution and feedback from stars is discussed 
in Sections [5] and [6l respectively. We then examine elliptical 
orbits of a gaseous stream in Section [T] and the question of 
the formation of mini star-clusters in Section [8] Finally, we 
summarize and conclude in Section [9l 



2 METHODS AND DISC FRAGMENTATION 
TESTS 

We use the SPH/iV-body code GADGET-2 
(ISpringel et al.. 200llSpringel. 20051 ) to simulate the 
dynamics of stars and gas in the (Newtonian) gravitational 
field of a point mass with Mbh = 3.5 x 10^ M©. The code 
solves for the gas hydrodynamics via the smoothed particle 
hydrodynamics (SPH) formalism. The hydrodynamic 
treatment of the gas includes adiabatic processes and 
artificial bulk viscosity to resolve shocks. The stars are 
mo delled as sink partic les, using the approach developed 
by I Springel et al. f 20051 ). modified in the ways described 
below. 

Table 1 lists some of the parameters and results of the 
simulations presented in this paper. The tests described in 
this section are those listed as S1-S5 (the "S" stands for 
"standard"). The units of length and mass used in the simu- 
lations are R\j = 1.2x lO^^cm ^ 0.04 pc, which is equal to l" 
at the 8.0 kpc distance t o the GG, and Mu = 3.5 x 10^ M©, 
the mass of Sgr A* fe.g.. lSchodel et al.. 20021 ). respectively. 
The dimensionless time unit is tu = l/r^(i?u), or about 

SPH particles are typically around 0.01 Solar masses (see 
Table 1). 



2.1 Gravitational collapse 

iToomre fl963 ) showed that a rotating disc is subject to 
gravitational instabilities when the Q-parameter 



2'kGp 



(1) 



becomes smaller than a critical value, which is close to unity. 
This form of the equation assumes hydrostatic equilibrium 
to relate the disc sound speed to its vertical thickness, H. 
Q = (GMbh/^^)^^^ is the Keplerian angular frequency, and 
p is the vertically averaged disc density. Gravitational col- 
lapse thus takes place when the gas density exceeds 



2nG 



(2) 



Ideally, numerical simulations should resolve the gravita- 
tional collapse down to stellar scales. In practice, this is im- 
possible for numerical reasons, and inst ead collisionle ss "sink 
particles" are commonly introduced ([Bate et al.. 1995) to 
model collapsing regions of very high density. To ensure that 
collapse is well under way when we introduce a sink particle, 
we require that the gas density exceeds 

Pcrit = PO + ^colPBH , (3) 

where po = 5 x 10~ g cm~ , and A 

col IS ^ large number (we 
tested values from a few to a few thousand). We in general 
find that our results are not sensitive to the exact values of 
po and Acoi, provided they are large enough (Section 15. 2p . 
Further details of our sink particle treatment are given in 
Section 13.11 



2.2 Gas cooling 

We account for radiative cooling via a simple approach with 
a locally constant cooling time. The cooling term in the en- 
ergy equation is 



© 0000 RAS, MNRAS 000, 000-000 



Star formation in a disc around Sgr 3 



\dt J cool 



tcoo\{R) ' ^ ^ 

where u is the internal energy of an SPH particle, and R 
is the radial location of the SPH particle, i.e. the distance 
to the SMBH. We parameterize the cooling time as a fixed 
fraction of the local dynamical time. 



tcool(R) — P X tdyn{R) 



(5) 



where tdyn = 1/^^ and /3 is a parameter of the simula- 
tions. This simple model allows for a convenient validation of 
the simulation methodology as such locally constant cooling 
time models were in vestigated in detail in previous litera- 
ture. iGammie f200lh has shown that self- gravitating discs 
are bound to collapse if the c ooling time, tcooi , is shorter 
than about 3/Q (i.e., (3 < 3). iRice et al. f 20051 ) presented 
a range of runs that tested this fragmentation criterion in 
detail, and found that, for the adiabatic index of 7 = 5/3 
(as used throughout this paper) , the disk fragments as long 
as < 6. 

Note that most of the simulations explored in this paper 
are performed for circular initial gas orbits, and also for 
relatively small disc (total gas) masses as compared with 
that of the SMBH. As will become clear later, this implies 
that during the simulations gas particles continue to follow 
nearly circular orbits, and thus their cooling time is constant 
around their orbits. One exception to this is a test done with 
eccentric initial gas orbits (Section [T]). 

2.3 Fragmentation tests 

To compare our numerical approach with known results, we 
ran several tests (S1-S5 in T able 1) in a setup reminiscent 
of that of iRice et al. f 20051 ). who modelled fragmenta- 
tion of proto-stellar discs. One key difference between 
marginally self- gravitating proto-stellar and AGN discs 
is the relative disc mass. Whereas the former become 
self- gravitating for a disc to central obj ect mass ratio of 
Mdisc/M* - 0.1 - 0.5 fe.g.. lRafikov. 20051 ). where M* is the 
mass of the central star, the latter become self- gravitating 
for a mass ratio as small as Mdisc/ Mhh — 0.003 — 0.01 (e.g., 
iGammie. 200llGoodman. 200^avakshin &: Quadra. 2003Levin. 
For marginally self- gravitating discs, the ratio of disc 
height to radius, H/R, is of ord er of the mass ratio, 
H/R ^ Mdisc/ Mbh toammie. 200lh . The disc viscous time 
is 



(6) 



whe re a ^ 0.1 — 1 is the effe ctive viscosity parame- 
ter (iLin k Pringle. IQSIG ammie. 2001) for self- gravitating 
discs. Thus, for H/R ^ 0.01, the disc viscous time is 
some 4-6 orders of magnitude longer than the disc dy- 
namical time. This is in fact one of the reasons why 
the discs become self- gravitating, as they are not able 
to heat u p quickly enough vi a viscous en ergy dissipation 
(|Shlosman Begelman. IQSdGammie. 200ll ). 

Also, since the disc viscous time is many orders of mag- 
nitude longer than the orbital time, as well as the total sim- 
ulation time, we expect no radial re-distribution of gas in 
the simulations, and indeed very little o ccurs. Our model s 
are essentially local, as emphasized by iNavakshin (20061 ), 
and hence it suffices to simulate a small radial region of the 
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Figure 1. Surface density of a region in a marginally stable disc 
from the simulation S4 (/3 = 4.5) at time t = 50. The map is cen- 
tred on the (x,y) = (1.1, —0.9) location. Many small and some- 
what dense gas clumps can be seen in the figure. None of the 
clumps however truly collapses, and most disappear on a dynam- 
ical time scale. 



disc. For definiteness, we simulate a gaseous disc of mass 
Md = 3 X 10^ M© ^ 0.01 Mbh extending from nn = 1 to 
^out = 4, where r is the dimensionless distance from Sgr A*. 

The total initial number of SPH particles used in runs 
S1-S5 is 4 X 10^. Gravitational softening is adaptive, with a 
minimum gravitational softening length of 3 x 10~^ for gas, 
and 0.001 for sink particles. The disk is in circular rotation 
nd a SMBH with mass of Mbh = 3.5 x 10^ M©. The 

height of H/R = 0.02 in 



disc is extended vertically to 
the initial conditions, which renders it stable to self-gravity. 
As the gas settles into hydrostatic equilibrium, it heats up 
due to compressional heating, and then cools according to 
equation Q. We used five different values of /3 = tcooi^^, in 
the runs, namely P = 0.3, 2, 3, 4.5 and 6 (see Table 1). 

The discs in runs with P = 0.3, 2 and 3 fragmented 
and formed star s, as expected based on the results of 
iRice et al. (20051 ). but our tests with P = 4.5 and 6 did 
not. However, the run with /3 = 4.5 and 6 did fragment in 
the sense of forming transient high density gas clumps, some 
of which can be noted in Figure [1] The maximum density in 
these clumps fluctuated between a few to a few dozen times 
Pbh, i.e. the clumps were dense enough to be self-bound, but 
not dense enough for our sink particle criterion, since Acoi 
was set to 30 for thes e tests (cf. Section [2T]) . 

iRice et al. f2005l ) reported fragmentation for values of 
P as large as 6. This difference in the results is most likely 
due to our discs being much thinner geometrically. The high 
density structures formed in our simulations, like filaments 
and clumps, are much smaller in relative terms than they 
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are in the simulations of proto-stellar discs. This is not sur- 
prising as the most unstable wavelength for t he Q ^ 1 disc i s 
of order of the disc vertical scale height H (jToomre. 1964 ). 
and our discs are much thinner. Also, the high density struc- 
tures are muc h more numerous in AGN discs. As analyti- 
cally shown bvl Levin f 2QQ6). self- gravitating clumps in such 
discs are quite likely to collide with each other on a dy- 
namical time scale. If clumps are tightly bound they may 
agglomerate. For tests with larger values of [3^ the clumps 
are marginally bound (as their maximum density exceeds 
pm by a factor of a few only) , and any interaction may un- 
bind them. In contrast, for high values oi H/R ^ Mdisc/M, 
as appropriate for proto-stellar discs, fewer clumps form, 
and the interactions are probably too rare to cause clump 
destruction. 



3 STAR FORMATION 

3.1 Two types of sink particles 

To design a physically reasonable and numerically sound 
procedure to deal with the gravitational collapse of cool- 
ing gas, we make use of well known results from the field 
of star formation. Simulations of collapsing gas haloes show 
that the first quasi-hydrostatic object to form has a char- 
acteristic size of a bout 5 AU, and may therefore be c alled 
a "first core" fe.g. jLarson. 196dMasunaga et al.. 1998 ^1. As 
this optically thick gas condensation slowly cools, the core 
decreases in size while accreting mass from the infalling en- 
velope. When the mass of the gas accumulated is of the or- 
der of 0.05-0.1 M0, H2 dissociation and hydrogen ionisation 
losses become efficient coolants, and the core can collapse 
dynamically to m uch higher densities and much smaller sizes 
([Larson, 196 6 Silk. 1977). The relevant time scale is of the or- 
der of ^10^ years, which is very short compared with free 
fall times of typical molecular clouds in a galaxy. Creation of 
the first collapsed objects is therefore nearly instantaneous 
in a "normal" star forming environment, and their sizes are 
sufficiently small so that they can be considered point masses 
compared with the cloud size of ^ 10^^ cm. 

The gas densities in self- gravitating AGN discs are 



p/mp 



10^ 



10 cm , i.e. many orders of magnitude 



larger than in galactic star forming regions, but smalleiQ 
than the gas density of a first core, ^ few xlO^^ cm~^. We 
assume that the initial stages of collapse should be simi- 
lar to "normal" star formation except that it starts from 
higher initial gas densities. For Sgr A*, at the distance where 
the young massive stars are located (Paumard et al., 2006'), 
the free-fall time in the collapsing gas clouds, tff ^ 
l/{Gpy^^ < tdyn, is of the order of 100-1000 years. For a 
more massive AGN, the time scales would be even shorter 
at comparable distances to the SMBH. The life-times of 
the first cores in this environment are then longer than the 
free-fall time, since the core collapse is dominated by the 
cooling time scale, which is much longer than the tff esti- 
mated above. This implies that these objects should actu- 



^ Closer in to SMBHs, the tidal shear density (pbh) becomes yet 
larger, and therefore there exists a minimum radius in the disc 
inside of which star formation does not take place, as proto-stars 
would be torn apart by tidal forces. 



ally be treated as having a finite life-time in our simula- 
tions. Furthermore, their sizes are only an order of magni- 
tude smaller than that of the disc scale height, H ^ 10^^ cm 
(|Navakshin Cuadra, 20051 ). Therefore, the fi rst cores hav e 
a decent chance to merge with one another ( Levin. 20061 ). 
They could also be disrupted by a more massive proto-star 
that passes by, form a disc and be subsequently accreted by 
the proto-star. 

In order to capture this complicated situation in a nu- 
merically practical way, we introduce two kinds of sink par- 
ticles. We call those with mass M < Mcore = 0.1 M© "first 
cores". Once the gas reaches the critical density (eq. [3]), we 
turn SPH particles into sink particles individually, so the 
initial mass of a first core particle is equal to one SPH par- 
ticle mass. First cores have geometrical sizes of i?core, and 
can merge with one another if they pass within a distance 
smaller than 2fmRcove of each other, where /m > 1 is a pa- 
rameter that mimics the effect of gravitational focusing in 
mergers. We assume for simplicity that the lifetime of the 
first cores is longer than the duration of the simulations. We 
performed a few tests with finite core life times and found 
very similar results; most of the merging of first cores hap- 
pens very quickly after they are created inside of a collapsing 
gas clump. 

Sink particles with mass M > Mcore are assumed to 
have collapsed to stellar densities, thus we refer to them 
simply as "stars". Stars are treated as point mass particles 
and thus cannot merge with one another, but they are al- 
lowed to accrete the first cores if their separation is less than 
Rcore- When the mass of a first core particle exceeds the crit- 
ical mass Mcore as a result of mergers or gas accretion, we 
turn it into a star particle. 

3.2 Accretion onto sink particles 

We use the Bondi-Hoyle formalism to calculate the accretion 
rate onto the sink particles. 



(7) 



(Ai;2 + ci)3/2 ' 

where M is the sink particle mass, and yo, Cs and are the 
ambient gas density, sound speed, and the relative velocity 
between the accretor and the gas, respectively. Once the 
accretion rate is computed, the actual gas particles that are 
swallowed by the sink particles are chosen from among the 
sink particle's neighb ors via the stochastic SPH method of 
ISpringel et al. f 20051 ). 

3.3 Maximum accretion rates 

As the gas density in these star forming discs is enormous 
compared to "normal" molecular clouds, the Bondi-Hoyle 
formula frequentl y yields super-Eddington accretion rates 
(Nayakshin, 2006*) . Assuming that gas liberates CM* / R^ of 
energy per unit mass, where M* and R^ are the mass and 
radius of the accretor, one estimates that the accretion lu- 
minosity is 



GM^M 



(8) 



where M is the accretion rate. Setting this equal to the 
Eddington luminosity, ivEdd = AnGM^mpc/ar = 1.3 x 
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1O^^(M*/M0) erg s ^, we obtain the Eddington accretion 
rate: 



(7T 

Note that this expression depends only on the size 
of the object, i?*. For first cores, i?* = i?core, 
which yields a very high accretion rate limit of al- 
most a Solar mass per year for i?core = 10^^ 
cm. For stars, we use the observational results of 
iDemircan Kahraman flQQllGorda Svechnikov flQQSh : 

/ \ 0.969 

— = 1.09 (^^j forM<1.52Mo, (10) 

/ \ 0.6035 

= 1-29 dnr) for M > 1.52 M©. (11) 

Hq \ Mo J 

For R = i?0, this yields an Eddington accretion rate limit 
of - 5 X 10"^ M© year"^ 

4 DISC EVOLUTION AFTER 
FRAGMENTATION 

We use the tests with cooling parameter (3 — 2 and 3, briefly 
described in Section 12. 3[ to point out some general trends 
seen in our simulations. These runs were continued far longer 
than was needed to form the first gravitationally bound 
clumps. Sink particles were introduced as described in Sec- 
tions [23] and [STT] The collapse parameter Aco\ was set to 30 
(see equation [3]) , and the gravitational focusing parameter 
is fm = 3. 

Figure [2] shows a column density map of a snapshot 
made at dimensionless time t = 75 for run S2 {[3 — 2). The 
left panel of the figure shows the whole disc, whereas the 
right panel zooms in on a smaller region centred on x — 1.8. 
Red asterisks show stars more massive that 3M0. We do 
not show lighter stars and first cores for clarity in the figure. 
Since gas dynamical and cooling times are shorter for smaller 
radii, high density regions develop faster at smaller radii. 

The sequence of events is as follows. The innermost edge 
of the disc forms collapsed haloes in which sink particles are 
introduced. These sink particles grow by mergers with other 
sink particles and gas accretion. With time, the gaseous disc 
is depleted while the stellar disc grows. 

Looking at larger R in Fig. [2l we observe an interme- 
diate epoch in star formation, where bound high density 
clumps have formed but there are no stars more massive 
than 3M0 yet. This is because some of these clumps have 
not yet satisfied the fragmentation criterion, and others sat- 
isfied it only recently, thus containing only first cores and 
low-mass stars. Finally, at the largest radii, we see dense fil- 
aments (spiral structure) but no well defined bound clumps. 

Figure [3] is identical to Fig. [2l but shows test S3 (i.e. 
[3 — instead of = 2). Comparing the two runs, it is 
apparent that fragmentation of the more gradually cooling 
disc (S3) is naturally much slower than in the test S2. There 
are fewer stars in the disc, and there are fewer high density 
clumps. As we shall see later, within the given model, this 
allows the stars to grow to larger masses, resulting in a more 
top-heavy mass function. 

Figure |4] shows the column density maps for the run 
S3 at two later times, t = 175 and t = 275, in the left and 



right panels, respectively. Star formation spreads to larger 
radii as time progresses. In the right panel of Fig. |4l the 
innermost region of the disc is almost completely cleared of 
gas by continuing accretion onto existing stars (rather than 
further disc fragmentation, see Section W72\ . At the end of 
the simulations that do fragment, i.e. S1-S3, very little gas 
remains. To characterise the rate at which star formation 
consumes the gas in the disc, we define the disc half-time, 
thaif , as the time from the start of the simulation to the point 
at which the gaseous mass in the disc is equal to one half of 
the initial mass. Table 1 lists thaif for each simulation. Longer 
cooling times allow the disc to survive longer. Nevertheless, 
even for the run S3, where thaif ^ ^dyn, the half-time of the 
disc is only around 12,000 years in physical units. 

One notices that some of the stars are followed and/or 
preceeded by low column depth regions. These are caused 
by angular mo mentum transfer from the st ars to the sur- 
rounding gas (jGoldreich &; Tremaine. 1980l ). If stars are 
massive eno ugh, they could open up gaps in accretion 
discs (|Sver et al , 199llNavakshin, 2004 ) in the way plan- 
ets open gaps in proto-planetary discs. We do not ob- 
serve such "clean" gaps around stars in our simulations. 
We believe the main reason for this is the large num- 
ber of stars that we have here. The interactions be- 
tween the gravitational potentials of neighboring stars de- 
stroy the resonant charachte r of the gas-star interactions 
(|Goldreich &: Tremaine, IQSOl V In addition, self-interactions 
between stars force stars on slightly elliptical and slightly in- 
clined orbits (Nav akshin Cuadra. 2005 ), which also makes 
gap opening harder. 



4.1 Evolution of disc vertical thickness 

The disc starts off as a razor-thin gaseous disc and evolves 
into a thicker stellar disc. This behaviour is to be expected. 
Initially, gas cooling is sufficiently strong, and the disc tem- 
perature is regulated to maintain a geometrical thickness of 
order of H/ R :^ Maisc/ Mph, which is equivalent to having 
Q ^ 1 (e.g., lGammie~ 001). When a good fraction of the 
gas is turned into stars, however, the rate at which the disc 
can cool decreases. At the same time, a stellar disc, con- 
sidered on its own, would thicken with time since stellar 
orbital energy is transferred into stellar random motions by 
the normal N-body relaxation processes. 

In the intermediate state, when both gas and stars are 
present in important fractions, disc evolution is not triv- 
ial. In the runs described in this section, stars transfer their 
random energy to the gas via Chandrasekhar's dynamical 
friction process (Navakshin, 2006). The disc swells as time 
progresses, as can be seen in the two edge-on views of the 
disc from simulation S2 shown in Fig. [5] . The gaseous disc 
can thus become thicker than it would have been on its 
own. However, depending on conditions (i.e. cooling param- 
eter [3^ initial total disc mass, mass spectrum of stars), the 
two discs can be coupled or decoupled. In the latter case 
the stellar disc has a larger geometrical thickness than the 
gaseous disc. In such a case the rate at which stars heat the 
disc is not trivially calculated. This is especially so if stellar 
radiation, winds and supernovae are taken into account, as 
stellar energy release above the disc is much less effective in 
disc heating than it is inside the disc. 
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Figure 2. Snapsliot of tlie disc column density at time t = 75 for test S2 {(3 = 2). Tlie gas is rotating clocl^wise in tliis and all the other 
figures in this paper. The left hand panel shows the full simulation domain, whereas the right hand one zooms in on a region of the disc 
centred dX x = 1.8. Stars with masses greater than SM© are plotted as red asterisks. 




Figure 3. Same as Fig.Obut for test S3 {(3 = 3), i.e. for a longer cooling time. Note the smaller number of stars and high density clumps 
in this test compared with Fig. [21 where the cooling time is shorter. 
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Figure 4. Snapshots of the disc for test S3 at times t = 175 (left panel) and t = 275 (right panel). Together with the left panel of Fig. [3] 
these snapshots trace the gradual depletion of the gaseous disc and the build-up of the stellar disc. 




Figure 5. Evolution of the gaseous and stellar disc thicknesses in the test S2. The image is centred at x = 1.3. The left panel is plotted 
at time t = 50, when ~ 20% of gas was turned into stars, whereas the right one is for t = 175, at which point only ~ 15% of gas remains 
in the disc. The disc thickens with time as the remaining gas is not able to provide sufficient damping and cooling of stellar A/'-body 
heating. 
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Figure 6. Disc fragmentation rate in the tests with j3 = 0.3 (SI, 
thick sohd), (3 = 2 (S2, dashed), and ^5 = 3 (S3, thin sohd). The 
shorter the cooHng time, the faster disc fragmentation proceeds. 

4.2 Fragmentation quenching 

We expect that, following an increase in disc effective tem- 
perature and geometrical thickness with time, the vertically 
averaged density of the disc must drop. The disc may then 
evolve into a non self- gravitating state as Q increases above 
unity (see equation [1]) . Further disc fragmentation should 
cease. This effect can be noticed in the right hand panel of 
Fig. [2] Near the inner edge of the disc, there are stars but 
no high density clumps or filaments implying that the disc 
is no longer fragmenting in that region. 

Figure [6] shows the radius-integrated fragmentation rate 
of the disc in the simulations with /3 = 0.3 (run SI, thick 
solid curve), (3^2 (S2, dashed) and /3 = 3 (S3, thin solid). 
The fragmentation rate is defined as the total mass of first 
cores created per unit time. Not surprisingly, the shorter the 
cooling time (smaller /3), the more rapid is disc fragmenta- 
tion. The more vigorous disc fragmentation explains why 
there are more stars and dense bound gas clumps in Fig. [2] 
than in Fig. [3] at the same time (t = 75). 

One can also see that the expectation of a decrease of 
the fragmentation rate with time is borne out. At the peaks 
of the respective curves, most of the disc mass is still in 
gaseous form. Therefore, the decline in the fragmentation 
rate with time is indeed mostly a consequence of a change 
in the disc state (higher Q-parameter) rather than due to 
the disc running out of gas. 

4.3 IMF of the formed stars 

It is interesting to compare the distribution functions of stel- 
lar masses from our simulations. The end state of our sim- 
ulations, i.e., when the majority of gas is turned into stars, 
corresponds to the "initial mass function" (IMF) of a stellar 
population. Therefore we use this name to refer to our mass 
distributions. Figure [71 shows the IMF of stars formed in 
the three simulations S1-S3. Table 1 lists the first two mo- 
ments of the distribution, i.e., the average mass, (M*), and 




0.1 1.0 10.0 100.0 

Figure 7. Mass function of stars ("IMF") formed in the simu- 
lations with (3 = 0.3 (SI, solid), (3 = 2 (S2, dotted), and (3 = ^ 
(S3, dashed). The dot-dashed power-law in this and the following 
plots corresponds to the Salpeter IMF. 

(^M^Y^'^ • is clear from both the figure and the table, that 
the longer the cooling time, the more top-heavy (or, equiv- 
alently, bottom-light) is the resulting IMF. This outcome is 
not surprising. As we saw in Section 14.21 fragmentation is 
fastest for the smallest values of the cooling parameter /3. 
Further, fragmentation stalls in our models when the stars 
heat up the disc above (3=1 (see Fig. [6]). At this point, disc 
fragmentation stops but accretion onto stars continues. The 
average mass of a star reached by the time the gas supply 
is exhausted is roughly inversely proportional to the num- 
ber of stars at the time the disc fragmentation stalls. As we 
find many more stars in the tests where fragmentation is 
quickest, these simulations will also have the least massive 
IMF. 

The shape of stellar IMF formed in our simu- 
lations is very different from the IMF in the Solar 
neighbou rhood, which is close to the power- law func- 
tion of ISalpeter floHBI l. Several authors argued that 
this might be expected based on theoretical gro unds 
(jMorris. 1993lNavakshin. 200dLarson. 200dLevin. 20061 ). 
The proportion of massive stars in the stellar population of 
these discs might be very important for the physics of these 
discs, e.g., their survival chances against star formation. 



5 SENSITIVITY OF RESULTS TO 
PARAMETER CHANGES 

Due to numerical limitations, our simulations do not include 
a number of important physical processes. In particular, ra- 
diation transfer and chemical processes taking place during 
collapse of a protostar are dealt with merely via prescrip- 
tions. These prescriptions employ parameters which we have 
tried to set to as physically realistic values as possible. It is 
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important to understand to which degree our results depend 
on the choices we made for the values of these parameters. 
Here we describe a selection of tests aimed towards this goal. 
These tests, labelled E1-E8 in Table 1, were done for slightly 
different parameters than runs S1-S5. In particular, the in- 
ner and outer radii were set to rin = 1 and rout = 2, respec- 
tively, and the gravitational focusing parameter was chosen 
to be /m = 1. The radiative cooling parameter f3 was set to 
3 for all of these tests. 



5.1 Maximum accretion rate 

We calculate the accretion rate using the Bondi-Hoyle for- 
malism, with the restriction that it is not allowed to exceed 
the Eddington limit (equations [71 and [9]). In principle, excess 
angular momentum in the gas flow around proto-stars could 
decrease the accretion rate below our estimates. To explore 
possible changes in the results, simulations El and E3 dif- 
fered only by the maximum accretion rate allowed, with 1 
and 0.1 times the Eddington limit, respectively. 

The resulting IMFs of these two runs are shown in 
Fig. [HI along with the IMF for test S3, which was already 
presented in Fig [T] Comparing the solid curve (simulation 
El) with the dashed curve (S3), we see that changes in rout, 
Aco\ and fm did influence the IMF, making it less top-heavy. 
The main effect stems from the reduction in the value of the 
gravitational focusing parameter, fm- Fewer mergers imply 
that the final average stellar mass is smaller. 

Comparing the runs El and E3, one notices that while 
the IMFs of the two tests are similar at the high mass end, 
around the peak in the curves, there is a larger discrepancy 
at the low mass end. The largest difference in the IMFs 
occurs at around the critical mass. Merit = 0.1 M©. The dif- 
ferences are nevertheless not as large as one could expect 
given an order of magnitude change in the maximum accre- 
tion rate. We interpret these results as follows. A star is first 
born as a low mass star in our simulations, and it then gains 
mass by either accretion or mergers with the first cores. If 
the maximum possible accretion rate onto the star is high, 
it "cleans out" its immediate surroundings of gas, and gains 
most of its mass by accretion. On the other hand, if the ac- 
cretion rate is capped at lower values, as in the run E3, then 
gas gravitationally captured by the star will first settle into 
a small scale disc around the star. This disc is somewhat 
smaller than the Hill's radius Rn = R{M^/?>MmY^^ ^ 
O.OIR for the star of mass M* = 10 M©. The density in 
that disc will often exceed the critical density for star for- 
mation, and hence new first cores will be introduced in close 
proximity to the star. These cores then mostly and rela- 
tively quickly merge with the star. Because the rate of gas 
capture from the larger disc (the Hill or the Bondi rate, e.g., 
[Goodman k Tan. 2004Navakshin k Cuadra. 20051 ) is con- 
trolled by the total mass inside the Hill radius, the rate 
at which the dominant star inside the Hill sphere is growing 
is comparable. The end result (the IMF for massive stars) 
is then similar, regardless of whether the captured gas was 
directly accreted or first turned into fluffy proto-stars (the 
first cores) and then accreted. 

Another difference between the results of the runs El 
and E3 is the longer disc life-time, thaif, for the former sim- 
ulation for rather obvious reasons. 




0.1 1.0 10.0 100.0 

M [M3 J 

Figure 8. IMF from the simulations with two different maxi- 
mum allowed accretion rates (see Section |5.1|) . i.e., 0.1 times the 
Eddington rate (E3 in Table 1; dotted curve) and the Eddington 
rate (El; solid). These curves are compared with the IMF from 
simulation S3 (dashed). 



5.2 Critical collapse density 

Our fragmentation criterium is controlled by two parame- 
ters, po and Aco\ (equation [3]) . The default values for these 
parameters are po = 5 x 10~^^ g cm""^, and Acoi = 60. 
These values were used in run E3. The simulations E2, E3a 
and E4 are completely analogous to the run E3 except for 
the two parameters controlling the fragmentation. In E2, 
the fragmentation density threshold is drastically relaxed, 
with po — 5 X 10~^^ g cm~^ and Aco\ — 9- In run E3a, 
= 5 X 10~^^ g cm~^, but Acoi = 60, as in E3. Finally, in 
the run E4 we used po = 5 x 10~^^ g cm~^, and Acoi = 300. 

Figure [9] shows the resulting mass function for the stars 
formed in these four tests. Evidently, a significant relaxation 
of the fragmentation criteria resulted in very large changes 
in the IMF at its the low mass end. The high mass end is 
not so strongly affected. In fact, the IMF of the tests E3 
and E4 are almost identical, implying that the results be- 
come independent of the value of the critical collapse density 
(equation [3]) when it is chosen large enough. 



5.3 First core size 

Runs E3 and E5 are different only in the size of the first 
cores, i^core (see Table 1). Figure [TOl demonstrates that the 
IMF of the simulations is remarkably similar at the high 
mass end, where it contains almost all of the stellar mass. At 
the low mass end, however, simulation E5 has significantly 
more mass than E3. Clearly, the smaller first core radius 
i^core permits fewer mergers and slower growth of first cores 
by accretion. The exact value of Rcore is hence important 
for the low-mass end of the IMF. 
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Figure 9. IMF of runs E2, E3, E3a and E4 (blue dash-dotted, red 
dashed, dotted and sohd curves, respectively). The simulations 
differ by the value of the critical collapse density (see Section 15.21 
for details). 



Figure 10. Comparison of the IMFs of runs E3 (solid) and E5 
(dotted). The simulations have a different "first core" radius pa- 
rameter, i^core = 10^^ and 3 X 10^^ cm, for E3 and E5, respec- 
tively. See text in Section [531 and Table 1 for details. 



5.4 SPH particle mass 

The simulations E6-E8 had identical parameter settings but 
used different initial SPH particle numbers to test the sen- 
sitivity of the results on the SPH mass resolution. Figure [TTI 
shows the IMF from these tests. Clearly, there is a significant 
difference between run E6 (solid; 1 million SPH particles) 
and the two others, E7 (dotted; 2 million) and E8 (dashed; 
4 million). There is a far smaller difference between the two 
runs with the highest resolution, though, suggesting that 
the results are reasonably close to full convergence. Given 
other uncertainties of the runs, e.g., the exact value of the 
parameter Rcore and the cooling parameter /3, etc., it seems 
sufficient to have an SPH mass resolution at the level of 
the run E7, which has msph = 0.02 M©. Note that except 
for the simulations E6 and E7, all the other runs presented 
here used msph = 0.01 M©, as in the test E8. 



6 IMF AND FEEDBACK 

iNavaks hin (20061 ) suggested that the luminosity of the 
young stars which accrete gas from within a star-forming 
disc is sufficient to heat up the disc, increasing the 
iToomre (1964 ) Q-parameter of the disc above unity and 
hence making t he disc stable to further fragmentation. 
INavakshin (20061 ) also suggested that this will lead to a 
top-heavy IMF for stars born in such marginally unstable 
self- gravitating discs. To test these ideas, we ran two addi- 
tional simulations that had thermal accretion feedback in- 
cluded. The feedback is implemented in the same way as in 
ISpringel et a l. (2005). with a parameter < F < 1 intro- 
duced to quantify the fraction of the accretion luminosity 
that is fed back as heat into the surrounding gas. The ac- 
cretion luminosity of course depends on the accretion rate 
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Figure 11. IMF of runs E6-E8. The simulations differ in the 
initial number of SPH particles used. In particular, A^sph = 1? 2 
and 4 x 10^ for E6, E7 and E8, respectively. 

onto the star (equation [8]) . The feedback luminosity, FLacc, 
is spread over the SPH neighbours of the accreting star ac- 
cording to their weight in the SPH kernel. 

Our two test simulations Fl and F2 had feedback pa- 
rameters of 0.01 and 0.5, respectively, and are otherwise 
identical to the run E3, which had F — 0. One notices right 
away from Table 1 that the disc lifetime thaif becomes longer 
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Figure 12. Mass function of stars ("IMF") formed in the sim- 
ulations without (sohd, run E3) and with feedback. The dashed 
curve is for run Fl, F = 0.01, whereas the dotted one is for F2, 
F = 0.5. Note that in simulation F2, only ~ 200 Mq of stars were 
formed by the end of the simulation. See text in Section [6] for 
details. 

for the tests with feedback compared with the test E3. In 
particular, even for Fl with F = 0.01, the disc lifetime is al- 
most twice larger than for E3. The test F2 in fact could not 
be ran until "completion" , i.e. until most of the gas turned 
into stars, as star formation was severely inhibited in that 
test. At the end of the simulation, after 10^ time units, the 
disc lost only about 10% of its mass. We hence could only 
estimate that thaif will be around 10^ time units. There- 
fore, the expectation of a stalling of star formation is con- 
firmed in our simulations with the large feedback parameter 
F. However, one has to note that in a more sophisticated 
approach the disc cooling time needs to be determined self- 
consistently, whereas we here kept it constant artificially. 

The IMFs of the runs Fl and F2 are compared to that 
of the run E3 in Figure [12] The IMF of the low feedback 
parameter case, Fl, is biased towards massive stars slightly 
more than that of the simulation E3, providing some sup- 
port for the predictions of lNavakshin (20061 ). However, sim- 
ulation F2 has not been run for long enough to produce a 
meaningful stellar mass function. In perspective, we suspect 
that the constant cooling time model is a poor approxima- 
tion especially when stellar feedback is explored. We shall 
develop a more realistic cooling model to test the influence 
of feedback on the disc evolution in the future. 



7 ELLIPTICAL ORBITS 

One interesting question raised by the observations of young 
stars in the central parsec of our Galaxy is that there is 
apparently a population of young stars on orbits that are 
quite eccentric (eccentricities of these stars are of the or- 



der of 0.7; see IPaumard et al.. 20061 ). The two-body relax- 
ation time scale is a s long as ^ 10^ years at this location 
(jFreitag et al.. 20061 ). so if the origin of these stars lies in 
in situ star formation, such large ec centricities could not b e 
built up from initial circular orbits ("Alexa nder et al.. 20061 ) ; 
consequently their initial orbits should have been eccentric. 

To examine whether stars can form in an eccentric 
gaseous disc or a stream, we set up an additional test that 
is identical to simulation El, except for the initial condi- 
tions. In this run, designated as "Ecc" in Table 1, the ini- 
tial gas configuration is that of a disc segment characterized 
by 4 < < 7 and azimuthal angle < < 7r/4. Fur- 
ther, instead of placing the gas on Keplerian circular orbits 
with V = vk{R) = ^{R)R, we decreased the gas velocity 
to V = 0.7vk{R) while keeping its direction perpendicular 
to its position vector and in the plane of the gas, as in a 
circular disc. These initial conditions are artificial, i.e. un- 
likely reached in a realistic infall of a gas cloud into the inner 
parsec of the Galaxy. However they allow us to test star for- 
mation in an eccentric gas disc under simple and controlled 
conditions. 

In addition, we also introduced in this simulation a fixed 
spherical potential of the olde r stellar cusp as deduced obser- 
vationally by Gen zel et al. (20 03). The total stellar mass of 
this cusp within the region being simulated is about ^ 10% 
of the mass of Sgr A*, i.e. small but not negligible. 

Figure [T3l presents several snapshots of the surface den- 
sity for test Ecc, along with stars more massive than 3 Mq. 
Initially the disc segment is torn into an eccentric spiral that 
makes a few rotations around Sgr A*. The spiral develops 
kinks (which would be essentially parts of spiral arms in 
a circular gas disc) that run at an angle to the gas spiral. 
Stars are born within these kinks, in orbits with eccentric- 
ities ~ 0.5, similar to t he original gas orbits. The process 
differs from the results of lSanders (19981 ) , who simulated the 
formation of young Galactic Centre stars on elliptical orbits 
in a sticky-particle approach. In that approach it was found 
that star formation occured in shocks due to self-collisions 
of the stream. It is possible that such a shock-mediated star 
formation would also occur in our model if gas was allowed 
to cool much faster than the local dynamical time. 

Due to the non-Keplerian potential in this test, gas 
and stellar orbits precess. The precession rate is different at 
different radii, and therefore gas orbits become mixed and 
somewhat circularized over time. This is most clearly seen 
in the last snapshot (lower right) in Fig. [TS] where an inner, 
only mildly eccentric gaseous ring is present. 

The simulation Ecc was run for twice longer than our 
usual 10^ time units. Nevertheless, at the end of the simu- 
lation only ^ 30% of gas was turned into stars. We hence 
estimate the disc half-life time thaif to be around 3000 time 
units. This is ^ 35 times longer than that of the correspond- 
ing simulation El with circular orbits. While some of the 
difference is simply due to a longer orbital time at the loca- 
tion of the gas in test Ecc, a fair fraction of the difference 
is due to a comparatively slower star formation. This is not 
altogether surprising, as gas is heated due to shocks in the 
eccentric, precessing disc of simulation Ecc. 

Figure [U] presents an edge-on view of the column 
density of the disc and stellar positions at the end of 
the run Ecc. It is quite noticeable that the stellar disc is 
much thicker than the gaseous disc. While it is difficult 
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Figure 13. Evolution of the eccentric disc simulation at times corresponding to t = 150, 450, 1000, and 2000, from left to right and 
top to bottom. Initially, the disc segment is tidally sheared into an eccentric spiral. While precessing, the spiral forms an eccentric disc. 
Outer and inner edges of the disc precess at different rates which leads to shocks and circularization of gas orbits. During the simulation, 
the outer ring of gas makes about 30 rotations and precesses around one full circle. 



to pinpoint the exact reason for this high geometrical 
t hicknes s, it is most likely due to interactions between star s 
(e.g., iNavakshin Cuadra. 2005lAlexander et al.. 2006h . 
which can be substantial given that this simulation has 
been run for over 10^ years. Disruptions of gas-star clumps 
at their orbits' pericentres or during collisions with each 
other in the disc could also be important. 



8 IN SITU FORMATION OF "MINI STAR 
CLUSTERS" 

Several authors speculated tha t self- gravitating AGN d iscs 
can form very massive objects. ICoodman Tan (2004l l ar- 
gued that massive stars formed in such discs will continue 
to accrete until they reach the "isolation" mass, Mi ^ 



-1/2 



Mi ^ 550 M, 









3/2 - 


104 Mo 





3 X 10^ Mo 



1/2 



(12) 



Such a super- massive star would presumably grow via 
accretion from a gas disc around it. The linear sizes of 
the disc would be of the order of the Hill radius of the 
star, Rh = R{Mi/3MBuy^^. The disc may become quite 
massive, and fragment further. In this way systems more 
complicated than a single massive star can be formed. It is 
not impossible that the compact star cluster IRS13E, which 

orbits S^r A* uncomfortably close, was formed in this way 

(|Milosavlievic Loeb, 2004Navakshin &; Cuadra, 2005lLevin, 

We searched for such massive bound systems containing 
many stars in the results of our simulations. While we did 
find many binaries of massive stars, nothing on the scale of 
the isolation mass was present. In terms of physical effects 
that could limit the growth of such groups, a too rapid disc 
fragmentation is the most likely one. Namely, for the smaller 
value of the cooling parameter (3, the disc fragments quickly 
into too many stars that then compete for the gas. No clear 
"winners" emerge, as can be seen from the shape of the 
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Figure 14. Edge-on column density and stellar positions for the 
last of the snapshots from Fig. [13] Note that the stellar disc is 
much thicker than the gaseous disc. 



IMF, which rolls over very quickly at the higher mass end 
(see Fig. [7|. 

For longer coohng times, i.e., /3 ^ 1, disc fragmentation 
is less vigorous, creating fewer stars, and leading to a more 
top-heavy IMF. In this case, while the gaseous mass is still 
comparable to the stellar mass in the simulations, we did find 
many binaries and multiple systems, some containing up to 6 
stars. However, closer to the end of these simulations, when 
the gas supply was exhausted, only tightly bound binaries 
remained. There seem to be two reasons for this. Firstly, 
while the system is still gas-rich, collisions between "mini 
star clusters" take place. The collisions tend to unbind the 
systems (since /3 > 1). A few examples of this could be seen 
directly in videos of the time evolution of the simulations 
produced from the snapshots. Secondly, when most of the 
gas is exhausted, individual multiple systems are likely to 
evaporate the less massive components in favour of tighten- 
ing the most massive binary. 

These results could in principle simply mean that we do 
not have a sufficient time resolution to integrate small scale 
dynamics of stars in these mini star clusters. To check this, 
we have ran an additional simulation, identical to the run 
S2, but in which the time step criterium was 2 times more 
stringent. In this simulation, not presented in the Table 1, 
individual particle time steps were on average half of what 
they were in S2. The result was very similar in terms of the 
IMF (i.e. the average stellar mass and the average stellar 
mass squared were different by ^ 3% and 10% only). How- 
ever, we did find a ~ 50% increase in the number of stars 
with several close neighboors. However, again, no group was 
even remotely close to the isolation mass scale. The largest 
group contained four neighboors. While this issue requires 
further numerical work, we believe that the absense of mas- 
sive densely packed stellar groups is not a numerical artefact. 



Nevertheless, independently of the precision with which 
the numerical integration is performed, the results might 
also strongly depend on the gas cooling physics. A more 
careful treatment of gas coohng physics is needed in order 
to test the model of in situ formation of IRS13E more thor- 
oughly. 



9 CONCLUSIONS 

In this paper, we presented some of our first attempts to sim- 
ulate star formation in an accretion disc around Sgr A* in the 
low disc mass case. Since the relevant physics is only we akly 
dependent on the mass of the SMBH ([Goodman, 2003l l, we 
expect that the results are relevant to more massive SMBHs 
as well. On the technical side, in our numerical method we 
used a locally constant cooling time prescription, and we al- 
lowed for a finite collapse time of stellar "first cores" and 
their mergers (Sections [2] and [3]). We limited the accretion 
rate of stars to a fraction of the Eddington limit. Our simu- 
lations verified that the code reproduces disc fragmentation 
in agreement with previous analytical and numerical work 
(Section 12. 3p . We also tested the dependence of the results 
on some of the plausible parameter choices (Section [5] and 
Table 1). 

The scientific results of our work can be summarized as 
follows: 

(1) Qualitatively, star formation in marginally self- 
gravitating discs, i.e. those that have cooling times just short 
enough to fragment, results in a top-heavy IMF. This is due 
to the suppression of fragmentation by disc heating either 
through N-body effects or feedback from the accretion lu- 
minosity. Mergers between first cores and gas accretion lead 
to a rapid growth of existing stars. Quantitatively, the IMF 
of our simulations does depend on physics that we do not 
completely resolve and which we model through several free 
parameters (see Table 1). 

(2) Without feedback, star formation is a very rapid process 
(unless cooling is inefficient, of course). As gas is bound to 
Sgr A*, nothing escapes from the star forming region, and 
essentially all of the gas in the disc is turned into stars. 

(3) Due to point 2, our models completely fail to feed 
the central AGN. Even the longest "living" gas disc in 
our runs (test F2 with feedback) would all be drained 
into young stars after about 10^ years, which is far 
too short compared with the normal timescale of feeding 
Sgr A* through viscous an gular momentum transfer (e.g., 
iNavakshin Cuadra. 20051 ). Modelling of higher mass discs 
is needed to establish whether star formation feedback in 
such discs could hold off the disc's demise due to star for- 
mation in favour of AGN accretion. 

(4) Discs more massive than marginally self- gravitating 
discs can be modelled in our approach with a small cool- 
ing time parameter, ^ 1. Without stellar feedback, such 
discs undergo a very rapid gravitational collapse. This leads 
to an IMF dominated by low mass stars. See point 5 below, 
however. 

(5) We have explored two limiting cases of marginally self- 
gravitating discs with thermal feedback. Putting even a 
small fraction of the stellar accretion luminosity back into 
the disc significantly increases the disc lifetime, that is, it 
slows down star formation. A more sophisticated numerical 
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treatment of heating/cooling is necessary to test whether 
the IMF of such discs will become top heavy, as predicted. 

(6) We also simulated star formation in an eccentric stream 
of gas in the gravitational potential of Sgr A* and the sur- 
rounding stellar cusp. The gas stream precesses, self-collides 
and forms an eccentric disc. Massive stars in eccentric or- 
bits are readily formed in the simulation. Tidal shearing and 
shocks in the eccentric case do not prevent collapse, though 
they do increase the lifetime of the gaseous disc. 

(7) None of the simulations formed an object similar to 
IRS13E, a bound "mini star cluster" in the GC. While more 
careful modelling is needed to confirm the result, it is sug- 
gestive of a formation mechanism for this object different 
from the in situ star formation in the disc. 

We thank Giuseppe Lodato for detailed comments on 
the paper. 
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Table 1. Description of simulations performed in the paper. Some of the run are printed several times in the table to allow an easier 
comparison between the tests grouped together. Notes: 

^ Cooling parameter introduced in Section [2?2l 

^ Feedback parameter, the fraction of stellar accretion energy fed back into the disc (Section [6]). 
^ Collapse over-density parameter (Section [JT} 

^ The maximum accretion rate of sink particles in units of their Eddington accretion rate (see Sections I5.1I &: [3?3]) . 
^ Gravitational focusing parameter for mergers of first cores (Section 13.11) 
Time scale for the gaseous mass in the disc to drop by half. 
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